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ABSTRACT 

The time delay between the formation of the progenitor systems of Type la supernovae (SNe 
la) and their detonation is a vital discriminant between the various progenitor scenarios that 
have been proposed for them. We use SDSS optical and GALEX ultraviolet observations of 
the early-type host galaxies of 21 nearby SNe la and quantify the presence or absence of any 
young stellar population to constrain the minimum time delay for each supernova. We find 
that early-type host galaxies lack 'prompt' SNe la with time delays of <100 Myr and that 
^70% SNe la have minimum time delays of 275 Myr - 1.25 Gyr, with a median of 650 Myr, 
while at least 20% SNe la have minimum time delays of at least 1 Gyr at 95% confidence and 
two of these four SNe la are are likely older than 2 Gyr. The distribution of minimum time 
delays observed matches most closely the expectation for the single-degenerate channel with 
a main sequence donor. Furthermore, we do not find any evidence that sub-luminous SNe la 
are associated with long time delays. 

Key words: supernovae: general; galaxies: elliptical and lenticular, cD; ultraviolet: galaxies 



1 INTRODUCTION 

The progenitor channels of Type la supernovae (SNe la) are of vi- 
tal importance both to our understanding of stellar evolution and 
\ to m odern cosmology where SNe la are u sed as standard can- 
dles dRiess et al.lll998r , |Perlmutter et alj|l999h . They are also a key 
part of our understanding of galaxy formation by virtue of their 
contribution to th e energy budget and chemical evolution of their 
host galaxies (e.g. | Gre ggio & Renzini 1983l: lMatteucci & Greggiol 

ll98d|Pipino&Matteuccill2004h . A range of recent observations 
have indirectly suggested that at least some SNe la can be produced 
by a variety of different progenitor systems, though these claims 
are tentative (Hansen 2003; Ruiz-Lapuente et al. 2004; Patat et al. 

l2007MVoss & Neleman s 2008; Ju stham et al.l2009h . Since no SN la 
progenitor system has been conclusively identified pre-explosion, 
we must use more indirect approaches to understand their origin. 

The various progenitor scenarios proposed in the literature 
have very different predictions for the time delay between the 
episode of star formation producing progenitor systems, and the 
time until the detonation of SNe la and so observational constraints 
on the distribution of time delays (DTD) are thus critical for dis- 
criminating between the various scenarios. We will discuss each of 
them in turn. 



1.1 Progenitor scenarios for SN la and their expected time 
delays 

Theory proposes that there are two main scenarios for the origin of 
SNe la; in both cases the system giving rise to the SN la is the ther- 
monuclear explosion of a Carbon-Oxygen (CO) white dwarf (WD) 
in a binary system. We briefly describe the main channels here; for 
a f ull review of progenitor sc enarios and explosion mechanisms, 
seelPodsiadlows kiet alj J2008h . 



1.1.1 Single degenerate scenarios 

In the single degenerate scenario (SD), the steady mass transfer 
from a main sequence (MS) donor star slowly builds the mass 
of a WD companion until it reaches approximately the Chan- 
drasekhar ma s p in it iating an SN la ( Nomoto & Kondo 1 199 it 
[ Hachisu et al.1 1 1999): [LangeretalJ |2000| : lHan & Podsiadlowskil 
boo4IJustham et all2 009). By the time of explosion, the donor can 
already be a slightly evolved subgiant. In this scenario, the time de- 
lay distribution is expected to pea k of ~ 670 Myr with virtua lly no 
contribution past ~ 1 — 1.5 Gyr ( lHan & Po dsiadlowski 200j). 

Both very short (or 'prompt') and ver y long ~ 1 Gyr time 
delay for SD channels have been proposed: [Hachisu et alj (2008) 



1 The point at which an explosive nuclear runaway occurs in a nomotat- 
ing CO WD may be slightly below the Chandrasekhar mass. N omoto et alj 
* E-mail: kevin.schawinski@yale.edu 1 1984) calculate a value of ~ 1.378 Mq. 
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suggest that mas sive (6-7 Mq) MS donors are viable, while 
IWangetal.l l 2009) argue for Helium star donors. In 

both cases, the detonation is expected to occur within < 100 
Myr. Very long SD time delays are possible if the donor is a low- 
mass red giant (RG) dHachisuetalll996l. Il999h . 

1.1.2 Double degenerate scenarios 

In the double-degenerate channel (DD) on the other hand, the ther- 
monuclear exp losion occurs when two white dwarfs in a binary 
system merge dlben&Tutukovlll984 ; Webbinkj [l984l) The range 
of possible time delays is set by the timescale for the formation 
of a WD binary system follo wed by the time needed for orbital 
decay via gravitational waves ( Shapi ro & Teukolskvlll983h . which 
depends only on the binary separation following the common en- 
velope phase. 

1.2 Early-type galaxies as an ideal laboratory 

Recent advances in our understanding of star formation - and ab- 
sence thereof - in the local early-type population, together with 
a large volume of new observations make early-type galaxies an 
ideal laboratory for constraining the time delay distribution of low 
redshift SNe la. Local early-type galaxies have generally been con- 
sidered to have formed at high redshift and devoid of any current o r 
recent star formation (e.g. lBower et al.iri992l ; lThomas et alj|2005l) . 
implying that any SN la progenitor systems are old with ages of 
several Gyrs. This view has been revised in recent ye ars by obser- 
vation s from the Galaxy E volution Explor er (GALEX; lMartin et al.l 
2005) UV space satellite: lYi et al. I « l2005t) found that a substantial 
fraction of early-type galaxies observed by GALEX had UV-optical 
colours incompatible with any scenarios of enhanced extreme hor- 
izontal branch stars (UV upturn; for a review of the current state 
see lYHl2008T) that could only be explained by the presence of small 
young stellar populations. 

This discovery was quantified using UV-optical colours to de- 
termine that between 1/3 and half of the local luminous early-type 
galaxy population harboured young stellar populations with mass 
fractions of a few percent and ages up to approx imately 1 Gyr 
dKavirai et alj2007] ; ISchawinski et al] [2006. 2007a) and their pres- 
ence has been detected out to z ~ 1 dKavirai et al. 2008). 

The inclusion of the GALEX near- and far-UV broadband pho- 
tomet ry (NUV, 1771-2831A; FUV, 1344-1786A; iMorrissev etal] 
120071) into the optical spectral energy distribution allows us to de- 
tect and quantify any young stellar populations in terms of age and 
mass-fraction, or to rule out their presence against the backdrop 
of an old, passively evolving bulk stellar population. This makes 
early-type galaxies in the local universe observed by GALEX an 
ideal laboratory to constrain and quantify the the ages of SN la 
progenitors for those events which occurred in early- type galaxies. 
Late-type galaxies with extended star formation histories are not 
suitable for constraining time delays out to interesting ages, as they 
always harbour stars of very young ages from ongoing star forma- 
tion. 

1.3 Previous constraints on SN la time delays 

Previous studies have endeavoured to determine the distribution 
of time delays (DTD) of SNe la and connect them to p rogenitor 
scenar ios. On the theory side, ana lytical calculations by iGreggiq 
d2005h and feelczynski etal] d2005h argue that the SD channel with 



short time delays i s unlikely to accoun t for all SNe la in passive 
systems. Similarly, |Pritchetetai]d2008l) make the case that the SD 
channel alone can only account for the observed relationship be- 
tween the SNe la rate and host g alaxy star formation rate under 
unrealistic as s umpti ons. Similarlv. lStrolger et al.l ( 2004, 2003) and 
iForster et al] d2006h attempted to re-construct the distribution of 
time delays from the cosmic star formation history, but they con- 
cluded that this approach is limited by uncertainties in the SFH. In 
addition, the small sample statistics o f lStrolger et ai] d2004l2005h 
are a further limiting factor. 

The properties of SNe la depend on the properties of their 
host galaxies. The brighter, slowly declining events tend to occur 
in star-forming host galaxies, while the dimmer, more rapidly de- 
clining SNe la are preferenti ally found in red, passive host galaxies 
(see e.g. lSullivan et alj2 006). This has led to the hypothesis that the 
bright SNe la result from progenitor systems with short time delays 
in star-forming galaxies, while the dimmer events, hosted by galax- 
ies without any current or recent star formation are due to scenarios 
with longer time delavs lMannucci et al] J2006h argued for the ex- 
istence of two separate SN la populations, a 'prompt' component 
with a time delay of ~ 1 00 Myr, and a 'delayed ' co mponent with 
time d elays of 3^1 Gyr. iBotticella et al] d2008l) and iTotani et"al1 
(2008) have analysed the host galaxies of SNe la and concluded 
that a substantial fraction of SNe la must have long time delays 
on the order of 2-3 G yr. At the other end of the delay time range, 
lAubourg et ai]|2008h have argued that there is a sub-population of 
SNe la with very short time delays of less than 180 Myr. 

iGallagher et al] ( 2008) recently argued to have found a direct 
link between the properties of individual SNe la and the bulk stellar 
populations of early-type host galaxies. They claim that that SNe 
la in host galaxies with older bulk ages and higher metallicities 
are dim mer and decay more rapidly. This has however been dis- 
puted by iHowell et al] d2009l) who find no such correlation. The 
tempting conclusion to draw here is that perhaps sub-luminous 
SNe la (sometimes called 1991bg-like, after the prototypical event) 
are from the long-time delay population, perhaps originating from 
a diff erent explosion mechanism due to the i r progenitor channe l 
(e.g. lHowellll200ll ; iTaubenberger et al]|2008l : IHowell et al]|2009l) . 
The consensus in the literature is that sub-l uminous events are 
due to explosions with low 56 Ni masses (e.g. lHdflich et~al]|2002l : 
iFesen et alj|2007l) . 

Rather than infer time delays indirectly from the bulk stellar 
population via optical data, in this Paper we probe deeper into the 
star formation history of a sample of 21 SN la host galaxies and 
use GALEX UV data to determine the presence or absence of minor 
episodes of recent star formation and thus constrain the minimum 
time delay for these SNe la. 

This paper is organised as follows. In Section [2] we discuss 
the selection of early-type SN la host galaxies that have been ob- 
served by GALEX and compare them to the general population. 
In Section [3] we motivate and describe our SED analysis method 
and present the results from applying it in Section [4] followed by 
a discussion in Section [5] We assume cosmological parameters 
(f2 m = 0.3, = 0.7, H = 70), consistent with the Wilkin- 
son Microwave Anisotropy Probe (WMAP) Third Year results a nd 
their combination of results with other data dSpergel et af]|2007h . 
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Figure 1. The UV-optical colour-magnitude relation of a sample o f normal early-type galaxi es, and the SN la host galaxies examined in this work. The grey 
points are the co mparison sample of earl y-type galaxies taken from Schawinski et al. (2007a), while the black points are the SN la host galaxies with GALEX 
observations from Forster & Schawinski ( 2008). For the SN la host galaxies, we show the error bars on the NUV — r colour; these are dominated by the error 
on NUV. We also label each host galaxy by the name of the supernova that occurred (c.f. Table[T}. 



2 SAMPLE SELECTION 

2.1 SN la early-type host galaxies 

We take the catalogue of visually inspect ed early-type galaxies 
that h ave hosted an SN la presented by iForster & Schawinskil 
(2008). The objects in this sample were culled by visually in- 
specting all host galaxies of SN la from the CfA list of super- 
novajjthat overlap with the Slo an Digital Sky Survey DR6 (SDSS; 
Adelman- McCarthv et alT2 008). The visual classification is vital to 
avoid the inclusion of systems with weak disks or spiral arms that 
can contaminate selections by proxies such as structural parameters 
or op tical colours (see e.g. Schawinski et al. 20 07bl : iLintott et al.l 
2008). We then limit this to those objects with observations in the 
GALEX GR4 archivtJE regardless of the survey; in total, this leaves 
us with a sample of 21 objects (see Table [T). We did not require a 
GALEX detection, only an observation, to avoid a bias towards UV- 
bright objects; however, all objects in our sample have a detection 
in at least one GALEX filter. 

Where GALEX data at multiple depths was available, the 
deeper image was used. The optical SDSS host galaxy and the cor- 
responding GALEX detection were matched by hand. 

2 http : / /www . cf a . harvard . edu/iau/list s/ Super novae . i 

3 http://galex.stsci.edu/GR4/ 




NUV-r 

Figure 2. The cumulative NUV — r colour distribution of both the SN la 
host galaxies (solid line) and the comparison sample (dashed line). The two 
colour distributions are closely matched. There appears to be a lack of very 
blue SN la hosts (NUV — r < 4), but neither a Kolmogorov-Smirnov 
test nor a Kuiper test indicate that the two populations are drawn from a 
different parent distribution; the lack of very blue SN la hosts is thus not 
statistically significant. 
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Table 1. The sample of SNe la and their host galaxy properties 



Galaxy Name 


Supernova 


Redshift 


NUV - r 


M r 


Best-fit 


Youn 


g age range 


Youngest possible t y 


Youngest possible t y 




Name 








young age t y 






at 68% confidence 


at 95% confidence 








mag 


mag 


Myr 


Myr 




Myr 


Myr 


NGC4493 


SN1994M 


0.02316 


5.69± 0.15 


-21.61± 0.002 


684 


650- 


1980 


719 


558 


CGCG169-00 


SN2002di 


0.03639 


5.57±0.13 


-21.98± 0.002 


881 


926- 


6670 


1193 


757 


IC4258 


SN2003an 


0.03704 


6.22± 0.13 


-22.36± 0.002 


. 


- 




6028 


2550 


NGC6095 


SN2003au 1 


0.03084 


5.96± 0.07 


-22.74± 0.002 


587 


531 - 


650 


558 


480 


NGC6109 


SN2003ia 


0.02954 


5.87±0.11 


-22.40± 0.002 


- 


- 




2424 


1025 


CGCG044-04 


SN2004bj 


0.05015 


5.68±0.17 


-23.20± 0.002 


757 


757 - 


7764 


926 


650 


NGC4493 


SN2004br 


0.02316 


5.69± 0.15 


-21.61± 0.002 


684 


650- 


1980 


719 


558 


CGCG089-05 


SN2004gs 


0.02664 


5.26± 0.04 


-21. 57± 0.002 


618 


531 - 


684 


558 


480 


IC0708 


SN2004H 


0.03168 


5.73± 0.05 


-22.64± 0.002 


- 


- 




4021 


1255 


2MASXJ1520 


SN20()5bm 


0.10350 


4.25±0.11 


-22.92± 0.003 


587 


433 - 


618 


505 


336 


MRK0693 


SN2005dh 


0.03836 


4.57± 0.03 


-22.98± 0.002 


618 


531 - 


837 


587 


480 


2MASXJ0141 


SN2005ex 


0.09000 


4.66± 0.13 


-21. 18± 0.006 


456 


336- 


531 


392 


275 


2MASXJ0134 


SN2005js 


0.07968 


6.19± 0.33 


-21.59± 0.004 


587 


433 - 


618 


505 


372 


2MASXJ01 10 


SN2005kt 


0.06540 


5.08± 0.10 


-21.00± 0.004 


684 


650- 


1389 


719 


587 


CGCG193-01 


SN2006bk 


0.04953 


5.89± 0.05 


-23.46± 0.002 


618 


558 - 


650 


587 


505 


NGC3841 


SN2006oq 


0.02120 


5.11±0.04 


-21.14± 0.002 


837 


719 - 


1320 


796 


684 


CGCG294-03 


SN2007ar 


0.05293 


6.06± 0.13 


-22.89± 0.002 








4449 


796 


NGC2577 


SN2007ax 2 


0.00686 


5.49± 0.04 


-20.30± 0.002 


757 


757 - 


6028 


837 


650 


CGCG077-10 


SN2007cf 3 


0.03293 


5.60± 0.07 


-21.41± 0.002 








1537 


975 


2MASXJ1410 


SN2007ei 


0.10000 


5.77± 0.46 


-22.04± 0.005 








1461 


837 


CGCG391-01 


SN2007jh 


0.04080 


7.08± 0.42 


-22.32± 0.002 








4021 


2191 



1 Spectrum resembles SN1991bg (Filippenko & Chornock 2003), potentially sub-luminous. 

2 Faintest SN la known fcasliwal et alj2008t). " 
A Reported sub-luminous bv lBlondin et alH2007l) . 



2.2 Are SN la host galaxies different from normal early-type 
galaxies 

Before we begin our quantitative analysis of the recent star forma- 
tion histories of the SN la host galaxies, we compare their UV- 
optical properties to those of normal early-type galaxies. For com- 
parison, we use the sam ple of early-type galaxies with GALEX ob- 
servations presented in Ischawinski et all ^20073) . This sample is 
magnitude-limited (r < 16.8), limited to 0.05 < z < 0.1 and vi- 
sually inspected for early-type morphology. We must preface this 
comparison by pointing out that while this comparison sample is 
well-defined, the SN la host galaxies sample is not; it is drawn 
from a heterogeneous parent sample of SN la composed of vari- 
ous surveys, all with their own selection effects and serendipitous 
discoveries of individual supernovae. Keeping this caveat in mind, 
we proceed to compare the two. 

In Figure [T] we show the NUV — r colour-magnitude dia- 
gram for the SN la host galaxies (black points) and the comparison 
sample from lSchawinski et al.l J2007a) (grey points). The compar- 
ison sample reasonably covers the parameter space of the SN la 
host galaxies. The host galaxies of SN2005dh and SN2005bm are 
both very blue in NUV — r and very luminous. We compare the 
NUV — r colour distributions of the host galaxies and the com- 
parison sample in Figure [2] While the cumulative distribution of 
the host galaxies might appear to indicate a lack of host galaxies 
at the bluest NUV — r colour, both a Kolmogorov-Smirnov and 
a Kuiper test indicate that the two distributions are consistent with 
being drawn from the same parent distribution, and so this lack of 
very blue SN la hosts is not statistically significant at the 95% level. 



3 METHOD 

3.1 Quantitative analysis of the recent star formation history 

The exquisite sensitivity of the near-UV to small amounts of young 
stellar populations makes it the ideal tool to quantify any young 
populations, or to rule out their presence down to very small lev- 
els. The old bulk stellar populations of massive early-type galaxies 
makes the exercise of measuring small young components easier. 
In Figure [3] (left), we show the evolution of the NUV — r colour 
as a function of the age of the young component for a range of pa- 
rameters to illustrate the strong age sensitivity of NUV — r. On the 
right, we show a composite SED of a single 12 Gyr old burst (red) 
and a 800 Myr young burst of a 1% mass fraction (blue). While 
it hardly changes the total SED (green) in the optical wavelengths 
probed by SDSS, the young population dominates the UV wave- 
lengths probed by GALEX. 

We parameterise the star formation history as two compo- 
nents: a burst on top of a old component. This parameterisation has 
been successfully impl emented before in order to answer questions 
of galaxy formation (Ferreras & Silk 200(J iKavirai et al.l 120071 : 
Schawinski et al. 2007b, 2009). Any young stellar population, even 
if small, is luminous compared to the underlying older population. 
Therefore, by marginalising over all possible star formation histo- 
ries before the most recent episode of star formation, we are able to 
constrain it, usually out to ages of ~ 1 Gyr and rule out any such 
episodes out to 1-2 Gyr. 

The properties of any young component is degenerate with 
a number of other parameters. Factors include the metallicity Z, 
which must be varied over the entire plausible range of metallic- 
ities, and should include an internal metallicity distribution with 
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Table 2. Model Parameters 



Parameter 



Range 



Old age 
Young age 1 
Metallicity 
Dust extinction 2 



to 
Z 

E{B - 



V) 



Young mass fraction f y 



1-15 Gyr 
0.1-15 Gyr 

i - 3.25Z Q 

0- 0.3 

1- 100% 



1 The you ng age t v is always re stricted to be less than the old age t . 

2 We use a lCalzetti etab 1 2000) extinction law. 



a tail ranging to v ery low metallicit i es. For old ages, both the 
low metallicity tail jPark & Leell 19971 : iMaraston & Thomasll200Cl) 
and hi gh metallicities can contribute moderate UV flux jYi et al.l 
1997al ltah. Since we do not wish to erroneously attribute any weak 
UV flux to young stars, we must account for this alternate, old ori- 
gin. The availability of far-UV data from the GALEX FUV filter 
aids somewhat in the differentiation, as these old populations tend 
to be hotter than young populations of a few hundred Myr and so 
have a steeper rising UV spectrum. We model our internal metallic- 
ity distri bution on the observed distribution of the bulges of nearby 
galaxies ISaraiedini & Jablonkj2005l : lHarris & Hamj200dl2002l ; 
Harris et all 1999). We range the metallicity Z such that the mass- 
weighted total metallicity ranges between ^Zp, — 3.25Zp 



20 © ~~ O.^OZjq. 

Further fa ctors are extinction by dust, which we model us- 
ing the law of Calz etti et al. I < l2000h and range the colour excess 
E(B — V) from 0-0.3 to account for levels of extinction typical 
in early-type galaxiefl and the age of the bulk old stellar popula- 
tion. We model this dominant component as a single stellar pop- 
ulation (SSP) of varying age from 1-15 Gyr. The young compo- 
nent is represented by a composite stellar population describing a 
brief burst of star formation declining with an e-folding time of 
100 Myr. The key parameters of this young burst are its age t y 
and it s mass-fraction /„ , which are somewhat degenerate with each 
other. ISchawinski et al.1 J2007bl) showed that the folding in of the 
luminosity-weighted age and metallicity measured from stellar ab- 
sorption indices (Lick indices) can break this degeneracy, but SDSS 
spectra are not available for all SN la host galaxies, and the major- 
ity of our host galaxies are so nearby that the 3" spectroscopic fibre 
of SDSS samples only the very central parfl which represents a bi- 
ased sub-sample of the bulk stellar population, as many early-type 
gala xies exhibit radial gradients in their ste l lar population para me- 
ters jDavies et alJl993l ; ICarollo et al.lll993l ; lTrager et ai1l200rj|) . 

The degeneracy between the mass-fraction and the age of any 
young stellar population represents a fundamental limit to the con- 
straints that can be derived from applying this method. The limit 
on the age of the young population t y that we can derive is a func- 
tion of the lowest mass-fraction f y that we probe down to. The age 
t y is not a strong function of the mass-fraction f y until below the 
~ 1% level, where even UV observations fail to be sufficiently sen- 
sitive. We thus range in mass-fraction f y from 100% down to 1%, 
far lower than is possible with optical photometry alone. All con- 



4 The measured range of E( B — V) for passively evo lving early-types is 
very low; see e.g. Figure 11 in lSchawinski et alj j2007bT) . 

5 In one case, the SDSS spectroscopic target selection algorithm selected 
two off-nucleus positions. 



straints on the age t y derived here are thus subject to the caveat that 
we do not account for young populations below the 1% level. 

We generate a library of photometric data points as a function 
of the parameters described above, resulting in 6.5 million model 
SEDs in the restframe of ea ch SN la host galaxy using the mod- 
els of lMarastor] d2005LI 1998b . We fit each host galaxy to all library 
SEDs and compute the \ 2 statistic for each, convert them to prob- 
abilities, and marginalise over all parameters except for the young 
age t y in order to obtain the probability distribution function Pity). 
In order to account for unknown systematic errors in the SDSS and 
GALEX zeropoints and other effects, we add in quadrature a uni- 
form error of 0.1 mag to each band. 

We have tested the effect of adding various amounts of sys- 
tematic error to broad-band photometry for similar codes in gen- 
eral, and for the data set and specific implementation used in this 
paper. By increasing the amount of assumed systematic error, the 
minimum set of parameters (t y , f y ) do not change and individual 
fits do not change from cases where a young population is detected 
to one where a young population is ruled out. Increasing the as- 
sumed error increases the size of the errors on (t y , f y ) at the min- 
imum, or lowers the limit on t y in cases where no young compo- 
nent is detected. While the SDSS photometric system is exquisitely 
calibrated, down to the 1-3% level for DRfJfl the photometric re- 
peatability o f GALEX is ~0 .03 an d 0.05 mag in NUV and FUV, 
respectively dMorrissev et alj|2007t) . Given this, we believe a 0.1 
mag error may be a conservative overestimate. We note that the 
main driver in detecting and ruling out the presence of young stel- 
lar populations are UV-optical colours whose observed range spans 
5-7 mag, substantially larger than the 0. 1 mag error. 



4 RESULTS 

We present the results of our SED fitting analysis in Figures|4]and|5] 
and TableQ] For each supernova and host galaxy in our sample, we 
plot the marginalised probability distribution function P(t y ) from 
100 Myr to 15 Gyr. The supernova host galaxies in our sample can 
be divided into two classes: those where a small young stellar popu- 
lations with ages of < 1 Gyr, are present, and those where any such 
young population is ruled out down to 1% by mass-fraction. This 
limit corresponds to the age resolution of the near-UV (c.f. Figure 
|3J. For those where a young component is resolved, we can calcu- 
late the best-fit age (Figure|4jl, while for those without, we can only 
report lower limits on the age of the youngest stellar population - 
and therefore SN la progenitor - present (Figure [5)- For compari- 
so n, the results for the general early-type population are presented 
in lKavirai etail j2007l) . 



4.1 SNe la host galaxies with detected young stellar 
populations 

For 14 out of 21 SN la host galaxies in our sample, we do detect 
and constrain the presence of a small (few percent) young stellar 
population in the host galaxy. The age of this young population 
represents the shortest possible time delay for these SNe la, though 
it is of course possible that these particular SNe have longer time 
delays, up to the age of the bulk stellar population. In Figure[4] we 



6 See http://www.sdss.org/dr6/ and Adelman-McCarth v et alj 
<20O8l) 
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Figure 3. The utility of GALEX UV data. Left: The evolution of NUV — r as a function of time elapsed since the last starburst on top of a 12 Gyr old 
passively evolving SSP. We show this curve for three different mass-fractions f y . All three models are extinguished by E(B — V) = 0.1. This demonstrates 
the enormous range in NUV — r colo ur for even smal l young populations and its strong dependence on the young age, and its weak dependence on the mass 
fraction. Right: A sum of a 12 Gyr old lMarastor] <2005l) single burst (red) and a 1% mass-fraction 800 Myr old CSP (blue) summing to a total SED (green). 
In the optical SDSS filters (indicated by their filter transmission curves), the addition of this small, young burst hardly changes the SED, while the difference 
becomes apparent and dominates in the UV filters probed by GALEX. 



show P(t y ) for all 14 host galaxies with resolved young compo- 
nentfl In each panel, we indicate the best-fit age with a dashed 
black line. We compute the 68% two-sided confidence interval and 
shade it with single gray hatches. In Table[T] we report the best-fit 
ages and the two-sided confidence interval. 



95% limit. The minimum age distribution ranges from 275 Myr to 
1 .25 Gyr, with two SNe at ~ 2 Gyr and has a mean age of 650 Myr. 
There are no cases where ages approaching ~ 100 Myr or below 
are allowed, so it appears that there are no 'prompt' SNe la in our 
sampl^E 



4.2 SNe la with no detected young stellar populations 

For the remaining 7 SN la host galaxies, we do not detect the pres- 
ence of any young stellar population down to the 1% level in mass- 
fraction. For these, it only makes sense to compute a one-sided con- 
fidence interval. We shade the 68% and 95% confidence limits on 
the P(t y ) distribution in Figure [5] and report these limits in Table 
[T] We also give these one-sided confidence levels for those host 
galaxies with detected young components. 

At 95% confidence, we find that four SNe la (SN2003an, 
SN2003ia, SN2004H and SN2007jh) have time delays of at least 
1 Gyr, with SN2003an and SN2007jh with minimum time delays 
in excess of 2 Gyr. Given the caveats discussed in Section [3] this 
implies that long time delays for SNe la do occur in nature and at 
least ~ 20% (4/21 ) of SNe la in early-type host galaxies have time 
delays longer than 1 Gyr. Due to the inclusion of GALEX obser- 
vations, the robustness of this statement is greater than that of any 
previous attempts to establish long time delays for some SNe la. 



5 DISCUSSION 

5.1 The distribution of minimum time delays for SN la 

We are now able to plot the distribution of minimum time delays 
(DmTD) for SNe la that occurred in early-type host galaxies. In 
Figure [6] we plot the histograms of the DmTD for SNe la at the 

7 We note that for three of these - SN2002di, SN2004bj and 2007ax - there 
is a substantial tail to ages of several Gyr and the best-fit age is around 1 
Gyr; these host galaxies straddle the limit of ages that are detectable. 



5.2 Sub-luminous SNe la? 

Several works have re marked that su b-luminous SNe la occur in 
early-type galaxies (e.g. lHowellll200lh . Under the assumption that 
these host galaxies contain only old stellar population, this was 
taken to be evidence that sub-lumino us SNe la have very long time 
delays on the order of several Gyr dHowellluOOll) . Although our 
sample of sub-luminous SNe la is very small, we find that in the 
case of two sub-luminous SN la, a shorter time delay is possi- 
ble. Our sample contains three sub-luminous SNe la, SN 2003au, 
SN2007ax and SN2007cf. O f these, SN2003au an d SN2007ax (the 
faintest SN la observed vet; lKasliwal et all 2008) both are in host 
galaxies with young stellar components, while only SN 2007cf 
has a 95% lower confidence limit of ~ 1 Gyr, i.e. only one sub- 
luminous SN la definitely has a long time delay. While this does not 
prove that sub-luminous SNe la have short time delays, it clearly 
leaves open the possibility that they do. 

5.3 Implications for progenitor channels 

What do our results imply for the various proposed progenitor 
channels? There are two ways we can explore implications of the 
results in this paper for progenitor channels. We can assume that 
the bulk of SNe la originate from the youngest stellar population 
present in the host galaxy, rather than from the old bulk population. 

8 The definition of what constitutes a 'prompt' time delay varies somewhat 
in the literature. While a value of ~ 100 Myr is common, others are also 
used. For example lSullivan et alj l200d) use 500 Myr. By their definition, a 
number of the SNe la studies here would be classified as prompt. 
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Figure 4. The probability distribution functions P(t y ) characterising the age of the most recent episode of star formation for those SN la host galaxies in our 
sample where a young component is resolved. We label each panel by the name of the supernova and its host galaxy. Note that SN1994M and SN2004br both 
occurred in NGC4493, so the P(t y ) is identical. The black dashed line indicates the best-fit age, while the single-hatched region shows the two-sided 68% 
confidence interval. 



This assumption implies that the DmTD in Figure[6]is a good ap- 
proximation of the true DTD and so we can take the DmTD and 
compare them to the predictions of various theoretical calculations 
of delay times. We stress that we cannot directly test with the cur- 
rent data whether this assumption is a good one and it remains 
plausible that the DmTD does not approximate the DTD well at 
all and that the young populations detected by GALEX are mostly 
unrelated to the observed SNe la. In the following discussion, we 



explore what the DmTD implies for progenitor channels of SNe la 
assuming that it closely approximates the DTD: 

5.3.1 Short time delay scenarios (< 100 Myr) 

Scenarios that give r ise to 'prompt' SNe la, such as massive MS or 
Helium star donors l lHachisu et al]|2008l : IWang et alj|2009h do not 
appear to occur in our sample. This is not necessarily obvious, as 
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Figure 5. The probability distribution functions P(t y ) characterising the age of the most recent episode of star formation for those SN la host galaxies in 
our sample where no young component is resolved. As in Figure|4] we label each panel by the name of the supernova and its host galaxy. The single-hatched 
region shows the one-sided 95% confidence limit, while the double-hatched region shows the 68% limit. 
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Figure 6. The distribution of minimum time delays for the SNe la in our 
sample. These limits are the one-sided confidence limits derived from the 
P(t y ) distributions in Figures H and H] (also given in Table [T}. The mean 
and median age are 800Myr and 650 Myr respectively. 



a few percent of early-type galaxies do have opti cally blue colours 
and host sufficiently young stellar populations (Schawins ki et al.l 



though the sample used in this Paper does not contain any 
such objects. 



5.3.2 Intermediate time delay scenarios (100 — 1000 Myr) 

The majority of SNe la in early-types have minimum time delays 
on the order of a few hundred Myr. If we assume that the actual 
DTD is approximated well by Figure[6] then the binary population 
synthesis predictions for the SD MS + CO WD scenario m atches 
the distribution of ages well (see lHan & Pod siadlowski [2004l) . 



5.3.3 Long time delay scenario (> 1 Gyr) 

Regardless of whether Figure[6]is a good representation of the ac- 
tual DTD, a minority of SNe la must have long time delays in ex- 
cess of 1 Gyr. The two progenitor channels that might account for 
these long time delays are the DP channel, can naturally pr oduce 
such long time delays dlben & Tutukov 1984; Webbink 1984)), and 
the SD channel with RG donors. Calculations for the rate due to 
WD+RG progenitor systems predict th at they cannot account for 
all SNe la in passive host g alaxies (e.g. lYungelson & Livid[l998l : 
lHan & Podsiadlowskiir2 004 ) 
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6 SUMMARY 

We have used SDSS optical and GALEX UV photometry to mea- 
sure the minimum time delay for a sample of 21 SNe la hosted in 
local early- type galaxies. We constrain the age, or rule out the pres- 
ence, of any young stellar population in the host galaxies down to 
the 1% level by mass. From this, we are able to construct the distri- 
bution of minimum time delays (DmTD) for the 21 SNe la in our 
sample. We find: 

(i) There are no prompt (< 100 Myr) SNe la in our sample. 

(ii) For 14 out of 21 SN la host galaxies, we detect and constrain 
young stellar populations, yielding a range of minimum time delays 
of 275 Myr to 1.25 Gyr, with a mean of 650 Myr. 

(iii) For 4 out of 21 SN la host galaxies, we can rule out the 
presence of any young stellar populations younger than 1 Gyr, im- 
plying minimum time delays longer than that. For two of these, the 
95% lower confidence limits rule out any time delays shorter than 
2 Gyr. 

(iv) Two out of three sub-luminous SNe la occurred in host 
galaxies with detected (few 100 Myr) young populations. Only one 
occurred in a host galaxy where young populations (< 1 Gyr) can 
be ruled out. 

We find that the predictions for the sin gle degenerate MS + CO 
WD channel faan & Podsi adlowski 2004) best matches the distri- 
bution of minimum time delays for the majority of our sample and 
that at least four require a scenario allowing for time delays greater 
than 1 Gyr, such as a RG + CO WD or DD channel. Under the ex- 
treme assumption that the DmTD is a good approximation of the 
true DTD, the SD MS + CO WD channel could account for up to 
70% of all SNe la in early-type galaxies, though the actual percent- 
age could be much lower in case this assumption is not warranted. 
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